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ABSTRACT 


Child and Langmuir introduced a solution to space 
charge limited emission in an infinite area planar diode. 


The solution follows from starting with Poisson’s equation, 








and requires solving a non-linear differential equation. 





This approach can also be applied to cylindrical and 


spherical geometries, but only for one-dimensional cases. 





By approaching the problem from Coulomb’s law and applying 
the effect of an assumed charge distribution, it is 


possible to solve for space charge limited emission without 





solving a non-linear differential equation, and to limit 





the emission area to two-dimensional geometries. Using a 


Mathcad worksheet to evaluate Coulomb’s law, it is possible 











to show correlation between the solution derived by Child 


and Langmuir and Coulomb’s law. 
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I. INTRODUCTION 


A. APPLICATIONS OF VACUUM DIODES 
Recent emphasis on missile defense and the expansion 


of electronic warfare has increased a need for high-energy 





electromagnetic weapons. Directed energy microwave weapons 
require high-energy electron beams to achieve the power 
levels adequate for standoff utilization, which can reach 
the gigawatt range for peak power, which is far beyond the 


power production of a semiconductor device. 





One method of creating electron beams is through 





vacuum diodes. In a vacuum diode, an electron cloud is 
created at a cathode and accelerated across a vacuum gap to 
a collector anode. This can be done through thermionic 


emission, field FFECE mitters, or explosive plasma 








emitters. For the pulse power requirements of high power 





microwaves, explosiv mitters ensure that enough electrons 
will be available to create currents necessary for high 


power devices within the desired turn-on time. 


Experimental work on diode emitters iS an ongoing 








subject, and is focused on diode geometries and materials. 
Suppression of plasma effects, cathode longevity, and the 


use of micro-fiber cathodes are all areas of interest in 





the creation of fieldable high-energy weapons. However, 





the electromagnetic effects of large current and charge 


densities in such a diode are not entirely understood. One 





effect is the screening process from the vacuum charges of 


the electric field produced by the diode. 





B. DEFINITION OF SPACE CHARGE-LIMITED EMISSION 


At high current densities, the electric field due to 





the electrons and/or plasma in the diode cannot be ignored. 


1 


Since the operation of these diodes is in the high current 


regime, this effect on the total electric field is 





important. 








In a parallel-plate vacuum diode, the electric field 


produced by the diode is a function of gap distance and the 





potential difference applied. The vacuum electric field is 








constant and uniform. In this case, the current would be 





solely limited by the capacity of the cathode to emit 





charge. However, xperiment shows that the current 





produced by a real diode is indeed finite, regardless of 





how much charge the cathode can emit. The reason for this 


is the decrease of the electric field due to the charge in 





the diode gap. 


As electrons enter the diode, the charge creates a 





field opposite to the field applied by the diode. Tf 





emission from the cathode is allowed to continue, the 








electric field of the charges can eventually cancel the 





diode field. At this point, the charges emitted by the 
cathode are no longer accelerated away from the cathode and 


current is limited. This is space charge-limited emission. 





The current produced by the diode is not limited by the 








ability of the cathode to produce charge, but by the 
canceling effects the current has on the total electric 


field in the diode. 


This limit has been explored since 1911. However, no 





closed analytic solution exists for finite diodes, which of 





course, all practical diodes are. Since ultimately, the 


power of an electromagnetic weapon is limited by the power 





produced by the electron beam, this limit is an important 


factor. 


II. CHILD-LANGMUIR LAW FOR SPACE CHARGE-LIMITED 
EMISSION 


A. GEOMETRY OF AN INFINITE VACUUM DIODE 
The derivation of space charge-limited emission is a 
corner stone of plasma physics. Space charge-limited 


emission is defined as driving the electric field at a 





cathode emitter to zero, preventing any increase in the 
current density. The first form of the well-known Child- 
Langmuir Law was published in 1911. Beginning with 
Poisson’s equation and the one-dimensional boundary 


conditions of a parallel infinite plate diode, the maximum 





current density can be found. 





The geometry of the parallel infinite plate diode 


consists of an infinite grounded cathode and an anode held 








at some potential, Vo. The two plates of the diode are some 
distance, D, apart. Current is then allowed to flow across 
the vacuum gap. When the space charge-limited condition is 





met at the emission surface, the electric field produced by 








the current in the gap will balanc th lectric field 











produced by the potential difference of the two plates. 





This geometry is reflected in the following diagram. 


@(D) =¥, 





J(z,r) 


(0) =0 








Figure 1. Diode Geometry 


B. DERIVATION OF CHILD—-LANGMUIR LAW 
The unknown quantity is the maximum current density, 


J(z,r) that will be produced by such a diode. A reasonable 





starting point is the Poisson equation, relating the 


Laplacian of the potential, @®, to the charge density, p. 


V@(z,r) = Le” (1a) 
é 


oO 





For th one-dimensional case, the potential does not 





vary in the radial directional, so the Laplacian reduces 





to: 


O’O(z) __ plz) 
Oz é 





(1b) 


oO 


The current density is also related to the charge 





density and velocity of the charge by definition. Since no 


charges are being created or destroyed in this steady-state 


problem, the relation is simply: 





J(Z) = p(Z)* v(Z) (2) 
Then the charge density is expressed as: 
J(z 
pz) = 22 (3) 
v(z) 


The conservation of energy is used to express the 
velocity in terms of known parameters of the diode: 
PE(0)+ KE(O) = PE(z) + KE(z) (4) 
Since the charge is assumed to be emitted with zero 


velocity, the kinetic energy at the cathode is zero. Since 





the cathode is grounded, the potential energy at the 


cathode is also zero, neglecting gravitation effects. 








Therefore, when expressions for the potential and kinetic 


energy of an electron are substituted: 
1 
O=-e* D(z) +m, *v(z)" (5) 


Solving for velocity yields: 


v(z) = [2 e@ (6) 


Substitution of the expressions for charge density and 





velocity into the Poisson equation results in an expression 








relating potential and current density. The expression is 


a second-order non-linear differential equation. 





O°O(z)__ Jz) [m, 1 
62" é 2e D(z) 


oO 





(7) 


The boundary conditions that lead to the space charge- 


limited solution are the potentials at the anode and the 








cathode, and the condition that the electric field at the 








cathode be zero. 


@(0) =0 (8a) 


@(D) =V, (8b) 


o0®(0) er, (8c) 











The last boundary condition is imposed because we ar 
solving for the space charge-limited emission case. Using 
the preceding boundary conditions, an outline of the 
solution to the differential equation follows. Multiplying 
both sides by the derivative of the potential and 





integrating: 
2 
(A _v/ ae + oo =", 1 (9) 
Oz OZ é, \ 2e az Jo 
The integrals are of the form: 


flevau)= [ua (10) 


Solving these integrals yields: 


2 1 
(2 _ wt Me #9*M24€ (11) 
Pa 6 é, \ 2e 


Here the constant of integration is zero from the 




















boundary condition that the electric field is zero at the 








cathode surface. Taking the square root of both sides: 


oo __, [J (m. oo ‘aay 
Oz &, \ 2€ 


Rearranging: 


1/4 
oan =2% [xf Me) xa (13) 
6, \2e 











Integrating this differential equation: 


1/4 
4 pt ign [Fa] Me) #24 (14) 
3 £.. deZe 


Again, the boundary condition of the grounded cathode 
forces the constant to zero. By applying the anode 
boundary condition, one can solve for J, the Child-Langmuir 


Space Charge Limit for emission density: 





A. (Oe yrs 
J.,=—-—€ *_o 15 
eae Wan D? oe 





This is the maximum current density that can be 





produce from a parallel infinite plate diode, given a 





potential difference, V., and spacing, D. Given this 
current density, substitution yields solutions to the 
following quantities: 


Potential: 


(z) =V (zZ) (16) 
. D 


Electric Field: 


1/3 
ei-Se%+(Z) (17) 








4 V z =213. 
gee oe 18 
Plz) 9° p? (=| (18) 


Cc. LIMITATIONS OF CHILD-LANGMUIR DERIVATION 
Although the Child-Langmuir approach solves for the 


space charge limit, ae is only valid for the one- 








dimensional case. It has been generalized to infinite 





cylindrical and spherical geometries, but it has not been 











solved for a finite emission area. Instead of approaching 
the problem from the Poisson eguation, which generates a 
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non-linear differential equation, if the integral form of 
Coulomb’s law is used, the space charge limit can be 


developed by integration. 


III. COULOMB’S LAW 


A. COMPARING POISSON’S EQUATION WITH COULOMB’S LAW 





In the diode configuration for the space charge 


limited formulation, it is assumed that the diode is in the 








steady state condition, and does not vary with time. With 


the addition of a background confining magnetic field, 





these conditions reduce Maxwell’s Equations for 
electrodynamics down to Gauss’s Law: 


VEG = 20 (19) 
é 


One result of Gauss’s law is the Poisson equation, 
from which the Child-Langmuir solution for space charge- 
limited flow is derived. As a definition, the electric 


field is the negative of the gradient of the scalar 











potential: 
E(r) =-VO(r) (20) 
Using the definition of the Laplacian operator: 
V -(-V®) =-V°@ (21) 
Substitution into Gauss’s Law gives the Poisson 
Equation: 
V*@(r) = ee) (22) 
é, 
Gauss’s Law also generates the integral form of 
Coulomb’s law. Substituting the integral of a Dirac delta 


function into the expression for p: 


pry 1 


ere [425° (r -r)p(r)dz' (23) 


One form of the delta function is: 


(rr) 























416° (r —r') = V -| ——— (24) 
(r=r') 
Substituting for the Dirac delta function: 
1 ~ = = lps] Gr) |_- 
— | 425° (r =r) p(r')dr' = —|V- LAP) | a@ydc (25) 
Ane Ane et 
¢ Q (r-r') 

Since the integral in the above equation is with 


respect to r', and the divergence is with respect to r, the 








divergence can be placed outside the integral, and the 
constant placed inside. From Gauss’s law, this is equal to 


the divergence of the electric field: 








(r-r) 


(r-r) 


5 le(r)dr' (26) 





VB) =V[— 








This must be hold for any r, so the arguments must be 








equal. Thus, Coulomb’s law is derived from Gauss’s law: 
= 1} c-h) | a. 
E(r) = | = \p(r)dr (27) 
47é - 
o| \(r—r') 








Therefore, Coulomb’s law and Poisson’s equation are 





two different, but equally valid forms of Gauss’s law. The 


results obtained by applying either equation to the 





conditions of the space charge-limited diode should be 
equivalent. 
B. COULOMB’S LAW GEOMETRY 

Coulomb’s Law in integral form is an extension of 


Coulomb’s force law for two point charges. Evaluating 











Coulomb’s law results in an expression for the electric 
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field produced by a distribution of charges. The integral 











form is: 
=-, lrpla,o 
E(r,) =——|— dq, ) (28) 
: ree ay Al 
In this case, th lectric field can vary depending on 
the position where it is measured, the field point, oe It 








is also dependent on the magnitude and position of the 


> 


charge distribution. The vector A is the displacement 


_ 


between the source point, r,, and the field point. 


A=(ry—r) (29) 


The geometry is reflected in the following diagram. 





dq 


Source Point 











Figure 2. Coulomb Diode Geometry 





The differential charge element, dq, can b xpressed 
as the product of the charge density, p, and a differential 


volume element. In cylindrical coordinates, this becomes: 








dq = p(r,)*r.dr,d0,dz, (30) 
To simplify the analysis of Coulomb’s law, the beam 
can be assumed to be confined horizontally by an external 
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magnetic field, thus require only the analysis of the 





vertical component of the electric field. In this case: 





Ae = (2) -%,)z (31) 


Although the directionality of ri and E have radial 








and azimuthal dependenc removed by the confinement 


assumption, it does not change the geometry of the problem. 





The magnitude of 4 remains: 





l= yr; +1, -2r,r,[sin(@, )sin(@,) + c0s(8, eos, |+ (2, 2, (32) 


Now that we have found expressions for all the 
components, the ability of software such as Mathcad to 


numerically solve integrals becomes attractive. 
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IV. MATHCAD WORKSHEET 





A. CHILD—LANGMUIR 

Mathcad is a commercially available numerical 
calculation program. Using Mathcad’s symbolic interface, 
it is possible set up a worksheet that compares the 


electric field predicted by Child-Langmuir and from the 








integral form of Coulomb’s law. If the two correspond, the 











flexibility of the worksheet will allow fOr the 


substitution for various charge distribution geometries. 





The first inputs into the worksheet are the global 





variables that control the initial geometry of the problem. 


For the initial purpose of comparing Child-Langmuir and 





Coulomb’s Law, the gap distance, Di; and potential 
difference, Vo, are constant factors, and can be set to 
unity to simplify the analysis. The same is true of the 
permittivity of free space, £&. 

The geometry of the Child-Langmuir diode is two 


infinite plates. The corresponding parameter iS an 





infinite radius, R. 





For comparison the results of the Child-Langmuir 


solution are defined. The functional form of the 





potential, @M(z), is given. Then the Child-Langmuir 











electric field, Ecl(z), is th one-dimensional divergence, 





or the derivative with respect to z. The vector component 





of the electric field in the z-direction is denoted by a 








subscript 2z. Ecl,(z) is plotted for comparison with the 


electric field calculated by Coulomb’s law. 
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—0.4 
—0.6 
Ecl 7(Z)_p g 
-1 
“12 
-1.4 
— 1.542 
0.2 0.4 0.6 0.8 
Z D 
Figure 3. Child-Langmuir Electric Field 
B. COULOMB’S LAW 


The charge density follows as the divergence of 
Ecl(z). This charge density will be used in the Coulomb’s 
law integral. However, a distinction must be made be in 


this transition between the source point, z;, and the field 





point, Zp In Child-Langmuir, the electric field is 





evaluated at the same z as the charge density. For 





Coulomb’s law, the charge density is integrated over all 
the sources to give the field at a distinct point. 


The Coulomb integral is integrated over all Oz, 





therefore the choice of 0: is arbitrary. Choosing 0: to be 


zero Simplifies th xpression for the displacement vector. 





The magnitude of the displacement vector becomes: 





lal = or, +7,7 2,1, eos@,)]+(z, - 2, ) (33) 


As a first check to the model, the conditions of the 
Child-Langmuir diode are applied. The symmetry of the 
infinite diode allows r¢ to go to zero. This: turther 
simplifies the displacement vector, removing azimuthal 


dependence. 





=. 2 2 
Mal = yr, +(z,-z,] (34) 
The evaluation of the Coulomb integral can be plotted 
as a function of Zg,. 
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0.667, 








—0.5 
P0067) Hy 
0 0.2 0.4 0.6 0.8 
0 Z f 1 
Figure 4. Coulomb Electric Field 
C:. COMPARING COULOMB’S LAW AND CHILD-LANGMUIR 
0.667, 
0.5 





Ecl 72) ' 
E A(z i 
E ee ‘) —0.5 


- 1.33324 5 


Z,Z¢.Z 5 D 


Figure 5. Comparison of Child-Langmuir and Coulomb’s Law 





It is evident from the plot that the limits of Coulomb 
field differ from the limits of the Child-Langmuir field. 





2 
Both the upper and lower limit are too high by -——. 





2V 
Subtracting 3D from the entire Coulomb field, it is seen 





that the Coulomb field differs from the Child-Langmuir 





field by this constant. 


The above plot shows that the Coulomb’s law integral 








and the Child-Langmuir solution differ only by a constant. 
D. ADJUSTING THE COULOMB FIELD 
That the Coulomb integral differs should not be a 





surprising result. The integral only accounts for the 





effect of the charge in the diode gap. However, the total 


electric field, as predicted by Child-Langmuir, is also 








influenced by the potential difference of the conducting 





plates. 

One approach to finding the effect of the conducting 
plates is the method of images. However, the solving for 
the images of a charge distribution on two conductors that 
are held at different potentials proved beyond the scope of 
this thesis. 

Another approach is to apply the boundary condition of 


the potential difference across the gap. 


D 
V,, = -[ | renee (zy )dz es 


Here, Exzcorrectea(Zt) 1S a sum of the electric field from 








the space charge, and a correction factor, the electric 


field produced by the conducting cathode and anode. 
| NEE (Z,) >= E, (Zy ) a5 jaan (Z,) (36) 











Integrating the space charge electric field: 
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D 1 
[, E.(z) dz, aera (37) 


Subtracting from the total potential difference: 


D 2 
i} | ere CaP )dz a ae (38) 


From this, the correction to electric field from the 





plates is: 


E ie (Zp) = - (39) 





This is exactly the magnitude of the difference of the 


Coulomb field from Child-Langmuir. This formulation shows 





a correspondenc between the predictions of Child- 
Langmuir’s solving of a differential equation and 
integration using Coulomb’s Law. 


E. GENERALIZING THE COULOMB INTEGRAL 





This integral formulation has the same limitations as 





the differential equation solution, in that it only holds 





for an infinite area diode. TO" b ffective for solving 








for a finite diode, the formulation must account for off- 


axis values of rez, introducing an azimuthal dependence. 





The generalized Coulomb integral does not have simple 


symmetries to allow for contraction of a dimension and is: 


(Z; =o) 





* 
ee re dr +r. ede i [cos(, )]+ (z, ye y 


2 
000 [ 1 +1, —2rr,Joos(6,)]+ &, -2,F } 








*r dé dr.dz, (40) 








In this form, Mathcad is unable to reliably solve for 








E,(Z¢) over the range O < Zz; < D. 
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V. CONCLUSIONS AND FUTURE WORK 


The goal of this thesis was to show a correspondence 





between the space charge limited conditions predicted from 


the Child-Langmuir law, and the conditions predicted from 





Coulomb’s law. In this respect, the Mathcad worksheet was 
successful. However, it should be noted, that assumptions 


imposed on the Coulomb integral limit the worksheet to the 








one-dimensional, infinite diode case. 


The assumptions that limit the Mathcad worksheet are 


assuming azimuthal symmetry of the charge distribution 





about the field point. For a finite distribution, this 





constrains the field point to be on the centerline of the 





distribution. This limitation does not allow for the 








exploration of th dg ffects of a finite beam, where the 








maximum current density exceeds the Child-Langmuir limit. 
These “wings” could vastly affect the performance of 
practical vacuum diode, as much more power would be produce 


than that predicted by the simple infinite approximation. 





Opportunity for future work exists in extending the 


worksheet to allow Mathcad to solve for arbitrary 





geometries. Including the azimuthal dependence would 
generalize the Coulomb integral, making the worksheet a 


powerful tool for exploring the electric field of a vacuum 





diode. 


9 





Hi 


IS PAGE 








NTENT 





ONALLY LEFT BLANK 





20 


APPENDIX 


The following pages are a rich text copy of the 


Mathcad worksheet used. 
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Global 
Subscripts ond and E represent vector 


Child-Langmuir: 


Potential: 
4 
3 


®O(z) = vo(7 
Electric Field: 


Ecl,(z) = d (2) 
dz i 


Ecl(z) > Sa 
Charge Density: 


p(,)= Poi Beles) =e 
P(Z, > 





Graph of the Electric Field: 


0.2 0.4 0.6 0.8 


Coulomb's Law: 
Infinite Diode Case: 
Displacement Vector: 





(% ~ zs) 
dol zp. 7%) = 5 5 
(r,) 7 (% = zs) 
(25, tit) = (r,)° + (2 - 24)” 
Charge Density: 
p(%) = = 
9-2. ? 
Electric Field: 
D -R 
il x 
way {| 
oO 
0 °0 


E,(.0003) = 0.581 





§ dependence removed due to 


azimuthal symmetry. 
Az, act) p(,) t dr, dz, 
(25, Zs ts) ) 
i 7 06 0.8 
ZF 
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E corrected zp) 5 E,(%) 7 : 





0.5 
Ecl,(Z) ; 
E. ze) 
E corrected 24)-0.5 
1 
-1.5 
: ue 0.4 ae 


Z,Z¢, Ze 


Space Charge Effect on Diode: 


= | ; Ea) Zp 


0101 


V = -0.338 
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0.8 


Allow for a Finite Diode: 


Displacement Vector: 


M 14.24.25 5585) = te + re - 2rp1,-(cos(0,)) + (2 - 25)” 


(4- 7.) 


Ia tp. 2p5t5 52595) = = s 5 
tee = 2r¢1,-(cos(0,)) + (2 - 7) 


z-component of the Electric Field: 





| (14.24% 075.06) 
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